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We analyze the leptonic sector in the left-right symmetric model dressed with a (Z2) 3 discrete 
symmetry which realizes, after weak spontaneous breaking, a small broken /1 o r symmetry that is 
suggested to explain observable neutrino oscillation data, r symmetry is broken at tree level 

in the effective neutrino mass matrix due to the mass difference m T 7^ m M in the diagonal Dirac 
mass terms, whereas all lepton mixings arise from a Majorana mass matrix. In the limit of a small 
breaking we determined #13, and the deviation from the maximal value of #23, in terms of the light 
neutrino hierarchy scale, m^, and a single free parameter h s of the model. 

PACS numbers: 14.60.Pq,12.60.-i,11.30.Fs 



I. INTRODUCTION 

The standard model (SM) has been a successful the- 
ory that explains most particle physics up to energies of 
about 100 GeV. However, as it is well known, it provides 
no explanation for the phenomenon of neutrino oscilla- 
tions that gives clear indications for the existence of non 
zero neutrino masses and mixings. Tiny squared neu- 
trino mass differences had been confirmed by a number of 
experiments, including Kamiokande, Super-Kamiokande, 
KamLAND, SNO, K2K, and Minos, among others. All 
known data are well understood in the setup where there 
are three mass eigenstates ^1,2,3, and three weak eigen- 
states ^ e ,^,T) which are related to each other through the 
mixing matrix Upmns PQ- This matrix, in the stan- 
dard parametrization, has three mixing angles named 
#23, #12, and #i3, one CP- violating Dirac phase, ip, and 
two Majorana phases. The mixing angles #23 and 812 
are measured in atmospheric and solar neutrino experi- 
ments, respectively, with very good accuracy, whereas the 
third angle, 6*13, and the CP phases have not been mea- 
sured yet, although the first is known to be small. Global 
data analysis on Ref. [5] provides sin 2 12 = 0.304+qqi6> 
sin 2 6» 23 = 0.50±8;g£, and sin 2 13 = 0.01±g;gjf, which is 
still consistent with zero ( see also Ref. [3] for an inde- 
pendent analysis regarding 6*13 ). A measurement of 613 
could be possible in future experiments [4] , though. Data 
also indicate that there are two scales for the squared 
mass differences, Am| = (7.65±%%l) x 10~ 5 eV 2 and 
Am 2 ATM = (2.4+°;i[ 2 ) x 10~ 3 eV 2 , that define the oscil- 
lation lengths at any given energy. All these prove that 
neutrinos are massive, and therefore that the SM needs 
to be extended to include their mass. 

The simplest route to include neutrino masses and mix- 
ings to the SM is to add the missing right-handed neu- 
trino (RHN) states to the matter content, and then in- 
voking the see-saw mechanism [5]. However, in this ap- 
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proach right the RHN mass scale is introduced by hand, 
with no relation whatsoever to the Higgs mechanism that 
gives mass to all other fields. The see-saw mechanism, on 
the other hand, comes in rather naturally in the context 
of left-right symmetric extensions of the SM [5J, aside 
from other nice features, as for instance the recovery of 
parity symmetry, and the appearance of right-handed 
currents at high energy, which makes such extensions 
also very appealing. Nonetheless, in any of these sim- 
ple matters or gauge extensions to the SM, a complete 
understanding of the peculiar mixing pattern of neutri- 
nos is not possible without additional assumptions. This 
is usually taken as a motivation to include flavor sym- 
metries to the models. As a matter of fact, a theoretical 
understanding of such numbers has been the goal of many 
theoretical works [7] in the last several years among which 
the bimaximal and tribimaximal mixing scenarios have 
attracted some special attention in the literature. Of par- 
ticular interest in some of those models is the fact that 
mixing angles are consistent with 6*23 ~ 45°, 6> 13 rj 0°, 
which, in the fcrmion basis where charged lepton masses 
are diagonal, seems to favor a fi r symmetry in the 
neutrino mass matrix. Indeed, a fundamental fj, •<-> r 
symmetry in the theory would exactly predict #23 = 45° , 
and 0i3 = 0°. 

Although fi o t symmetry is indeed a flavor symmetry 
in SM interactions, it is not an exact one for the whole 
particle physics due to the different masses of muon and 
tau leptons. Nevertheless, having two of the predicted 
mixing angles so close to the observed values could be a 
good indication that ji •<-> t symmetry does have some- 
thing to do with fixing the observed lepton mixings at 
some level. As it was shown in a recent analysis [8], the 
effect of the breaking produced by the sole charged lep- 
ton sector on neutrino mixings is rather negligible. This 
is because in a SM with three right-handed neutrinos, 
dressed with an effective [i O r symmetry in all neu- 
trino couplings, the \i and r mass difference enters into 
the mixings only through one-loop quantum corrections 
mediated by the W boson, and thus, it comes out very 
suppressed. Besides, the implementation of the fi -H> t 
symmetry in the SM turns out to be quite subtle, since 
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RHN's do not lie along similar representations as other 
leptons. RHN's are rather singlets of the theory, which 
drives us into an ambiguous definition of the meaning of 
lepton flavor in this sector. As a consequence, there is 
no unique or natural way to associate the properties of 
RHN under fi •<-> r symmetry. Only two generic classes 
of models do actually exist, depending on whether or not 
RHN's carry lepton number with them, though. 

The situation is expected to be different in the left- 
right symmetric model (LRSM), based on SU(2)l x 
SU(2)r x U(1)b-l gauge symmetry, with explicit parity 
symmetry LRSM has some advantages over the SM, in 
particular, as we already commented, it includes auto- 
matically right-handed neutrino fields in its matter con- 
tent, providing a working frame where left-right handed 
fields are treated on the same footing. This is partic- 
ularly important to give meaning to the mass scale of 
the right-handed neutrinos, which here is now linked to 
the breaking of parity, as well as to the scale of right- 
handed weak interactions. LRSM is also known to allow 
for Dirac and Majorana mass terms for neutrinos, and so, 
the see-saw mechanism naturally arises from it. A com- 
mon draw back of the model is that it introduces many 
additional free parameters which have to be constrained. 
However currently this is hard to do due to the lack of 
experimental data. All the phenomenology of LRSM is 
so far only constrained by SM data. An extended analy- 
sis of the LRSM parameters has been done in [5] . Flavor 
symmetries in the context of LRSM had also been previ- 
ously explored (see for instance 10J) and, as is usual on 
flavor extensions, extra particles are needed to support 
such symmetries. 

In the minimal LRSM Dirac neutrino masses would 
now be subjected to obey similar patterns than those 
on the charge lepton sector, due to the common origin 
of mass terms, which usually arise from the same set of 
Yukawa couplings. This is of course a consequence of 
the fact that in the LRSM, RHN's are paired with the 
right-handed leptons, which introduces by construction 
a meaningful lepton number for both the handed sec- 
tors. This makes the identification of/if>T symmetry 
neat. The realization of flavor symmetries becomes more 
constrained and a bit more challenging from the theo- 
retical point of view, though. In the particular case of 
fi •<-> t symmetry, its realization requires the model to 
provide diagonal charge lepton masses at tree level, from 
where different /i and r masses should arise naturally. 
This comes with the consequence of having a similar (di- 
agonal) mass structure for the tree level Dirac neutrino 
masses, and thus, introducing an extra tree level source 
for the breaking of fi r symmetry. Yet, one would 
like to have small predictions for the value of #13 and the 
deviation of maximality for 9 2 3, at most closer to the cur- 
rent experimental limits. The smallness of the induced 
breaking terms would possibly be a matter of parame- 
ter tuning, though. Another generic implication of these 
models would be that Dirac masses should not generate 
mixings. Lepton mixings would come solely from the Ma- 



jorana mass terms, reducing in this way the number of 
free parameters that the theory would involve to repro- 
duce the low energy data. To explore these aspects of 
such models is the main task of the present work. 

/i f-> r flavor symmetry can be implemented in any 
model with the help of discrete groups like Z2, S3, and 
A 4 among others [IT], which still makes it possible to 
build different realizations of the model with exact and 
broken /1 o r symmetry . We will follow this route for 
our model with the goal of providing a framework where 
a small broken fi <-> r symmetry is realized and where 
charged lepton mass matrix wig is given in a properly di- 
agonal form at tree level. Thus, we propose a LRSM with 
a (Z 2 ) 3 set of discrete flavor symmetries that provides the 
realization of a slightly broken fj, t symmetry, bro- 
ken only by the diagonal Dirac mass terms at tree level, 
with lepton mixings arising from the (yet fi ■<->■ r sym- 
metric) Majorana mass matrix alone. ( For models that 
use similar discrete symmetries in the context of the SM 
see, for instance Ref, [12].) In this model, the breaking 
of the discrete symmetries shall induce sizable non zero 
#13, and the deviation of O23 from 45°, controlled by a sin- 
gle parameter associated with the Yukawa couplings on 
charged lepton and Dirac neutrino masses. The paper is 
organized as follows: first of all, we present our extended 
LRSM with discrete symmetries and briefly argue on the 
scalar sector in Sec. II. In Sec. Ill, we comment on the 
way we realize a small broken fi «-> t symmetry. We 
present, in Sec. IV, our predictions for the mixings with 
a minimal number of free parameters, assuming no CP 
violation. Finally, we show the results obtained in the 
above section, in the presence of CP-violating phases, do 
not change drastically in general. We close our discussion 
with a summary of conclusions. 

II. THE EXTENDED LEFT-RIGHT MODEL 
WITH h^t SYMMETRY 

The model we shall consider along our discussion is 
based on the usual SU(2) L x SU(2) R x U(1) B -l gauge 
symmetry, but with an enlarged matter content that in- 
cludes two scalar bidoublets 4>i(B — L = 0) (i = 1,2), 
and four SU(2) R triplets, Aj (B — L = -2) (j = 1,...,4), 
aside from the usual left L a (B — L = —1) and right 
R a (B — L = 1) leptonic doublet representations, for 
a = e, /i, and r. Here we assume parity conserva- 
tion, where parity means the transformations: L R, 
4>i f-> 4>\. Thus, there would also be a similar set of 
SU(2)l triplets. Nonetheless, they would not be involved 
in the mass generation since they would have zero vac- 
uum expectation values (VEVs), and thus, we will not 
write them down explicitly. 

In order to get a neutrino sector that is compatible 
with observations we will introduce three Z 2 discrete 
symmetries under which the matter fields should trans- 
form as: 

Z 2 : (j) 2 ^ -h', L ^ <"> L r\ <-> Rt\ A 2 <-> A 3 ; 
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Z' 2 

z'i 



; L„ <-> -L T ;Ru <-> -R T ; A 2 O -A 3 



L u — L u ; i? u O — ii u ; A3 



-A,; A. 



where other fields not indicated remain invariant. Notice 
that Z2 realizes fj, •<-> r symmetry for the model in an 
effective way. It mimics \i •<-> r symmetry for the lep- 
ton sector, provided that ^2 is an odd field under such a 
symmetry. In the broken phase, however, the model will 
provide a mass pattern that will certainly correspond to 
a slightly broken fi ■<->■ r symmetry, as we will see. In ad- 
dition, Z 2 and Z'2 are additional symmetries whose role 
is to cancel all off-diagonal couplings among charged lep- 
tons, and thus the appearance of a flavor-changing neu- 
tral current at tree level. The above class of symmetries 
has been used before in the context of SU(2)l x U(1)y 
models with three vr right-handed neutrinos |12) . but in 
such a model fx ■<-> r symmetry remains exact, and lep- 
ton mixings emerge from the Majorana mass matrix. In 
our case, we will rather obtain a framework where this 
symmetry is broken at tree level in the effective neutrino 
mass matrix, following the spontaneous breaking of sym- 
metries. Some details regarding the analysis of the scalar 
potential with this particle content and symmetries are 
given in the Appendix. 

Next, according to the above symmetries, the most 
general Yukawa couplings of leptons to Higgs bidoublets 
are given as 



h(f>i + h(j)\ 



Re + L h 



ft-i^i + /ii0i — h 2 (j)2 
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R„ 



L T h\4>i + hi4>i + h 2 4>2 + h 2 4>2 Rt + h.c 



where <j) — t 2 0*t 2 . Parity symmetry is explicitly shown 
in the above equation. This is a feature of the LRSM 
before spontaneous symmetry breaking. On the other 
hand, we can see that the combined discrete symmetries 
restrict enough above couplings, in a way such that only 
diagonal mass terms arise. Then, we can identify our 
basis with the true flavor. As is clear, after symmetry 
breaking one shall get the correct charged lepton masses 
7^ m T , signaling the breaking of fx -H- t exchange 
invariance. By construction, such a breaking will also 
appear in the Dirac neutrino mass terms, which shall 
follow a similar profile. Moreover, this will be the only 
place where the exchange symmetry is broken afterward. 

The Lagrangian for the Yukawa couplings to triplets, 
allowed for the symmetries, is written as 



h{R e ) c R e + h {(R») c Rn + (Rt) c R t 
+/12 ' _ 
+A 2 R T 



A x 



((^) C A 3 + (R r ) c A 2 ) R e + (R e ) c ( A 3J R M 

- /23 



(R^) C A 4 R T + (R T )CA 4 R^ 
[R^> L] + h.c. 



(1) 



Here, the four triplets help us build a Majorana matrix 
which generates lepton mixings. As a matter of fact, since 



the Dirac mass matrix will be diagonal, the only source of 
mixings should arise from Eq.([I]). It is worth noticing the 
way the model separates the physical sources of mixings 
from which breaks //Or, although both components get 
finally interconnected by the see-saw at low energy. 

To get masses for all fermions, we have to break the 
gauge SU(2) L x SU(2) R x U(1)b-l symmetry in the 
usual way. So, the scalar field vacuum expectation values 
are given as 



h 
k' 



and 





Vj 



(2) 



where, in order to have the see-saw mechanism [5] at 
work, we have to assume that right-handed triplet VEVs 
are larger than those of the bidoublets. So far, for sim- 
plicity we are leaving out CP-violating phases in our anal- 
ysis, and so we neglected all phases in the vacuum ex- 
pectation values, as well as in Yukawa couplings. We do 
remark that this assumption is motivated by the larger 
number of scalar fields in the model, and it is taken in 
order to reduce the free parameters of the model to be 
considered, in this case CP-violating phases. Neverthe- 
less, in general, we do not expect CP phases to drastically 
change the general predictions of the model. Particu- 
larly, we see no indication that the order of magnitude 
of sin# 13 and of 1/2 — sin 2 623 could be affected in a sig- 
nificant manner, so we rather prefer to keep the analysis 
simple, since estimating such numbers is the main goal 
for what follows. In addition, to give enough support to 
our statement we will explore, later on, the model in the 
presence of CP-violating phases in an effective way. 

As we already mentioned, the charged lepton mass ma- 
trix, mi, is diagonal as well as the Dirac neutrino mass 
matrix, ran- Their entries are explicitly given as 

m e = hk-^ + hki, m^^ = h 1 k 1 + hiki^ (h 2 k 2 + h 2 k2); 
fh e = hki + hki, rh^ jT = h\k\ + h\k 1 =p (h 2 k2 + ^2^2) 

(3) 

respectively, where Dirac neutrino masses, wig , have been 
labeled according to the corresponding lepton number. 
Notice that the model produces a /i-r mass difference in 
a natural way. Notice also the similar profile that arises 
in neutrino masses where ll-t symmetry turns out to be 
broken by the essentially the same source. An interesting 
point that arises here is that the hierarchy between the 
last (m M and m T ) two Dirac masses is much softer than 
the {rriu and m T ) charged lepton masses, as a result of 
Eq. pj. Even if there were strong cancellations among 
the first and second terms in to m , as required to get the 
right hierarchy of charged lepton masses, that would not 
be the same for m M , which allows for some level of de- 
generacy on Dirac neutrino masses, in accordance with 
fj, •<-> r symmetry. At this level parity symmetry is bro- 
ken since we have taken zero VEVs for the left-handed 
scalar triplets, and therefore, the type-I see-saw mecha- 
nism arises to explain small neutrinos masses. 
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On the other hand, we do stress that the right-handed 
neutrino mass matrix, Mr, is symmetric under /i O r. 
This is due to the fact that the vacuum expectation 
values that minimize the potential satisfy A2 = A3 
(V 2 = V3). This happens as a consequence of the dis- 
crete symmetries of the model, as we explicitly show in 
the Appendix. Thus, we get 



AVl f 12 V 2 f 12 V 2 

Mr = I f 12 V 2 f 2 Vl /23V4 

fi 2 V 2 f 23 V 4 f 2 Vx 



(4) 



III. MASSES AND MIXINGS 

The matrix M v is diagonalized by Urmns ( m the stan- 
dard parametrization) , with 023 ~ a — f where \a\ <C 1. 
Notice that the a parameter can lay in the first or fourth 
quadrant of parameter space. 9i 2 is then identified as 8q 
and now #13 should be non zero. When 6* 13 and a go 
to zero, the mixing matrix becomes that of the bimaxi- 
mal framework. Under this approximation, the following 
three conditions should be fulfilled in order to get a di- 
agonal mass matrix: 



which clearly exhibits a fi-r symmetric structure. It is 
worth mentioning that even though Mr is the only source 
for neutrino mixings, by itself alone it does not fix #23, 
nor ^13. It is the conspiracy with the Dirac mass terms 
that do the fixing. For simplicity, we parametrize the 
inverse Majorana mass matrix as 



Ml 




(5) 



which of course maintains a symmetric profile. These 
new parameters are easily rewritten in terms of the Mr 
elements, but for the low energy analysis it is enough to 
use them instead, without obscuring the physical mean- 
ing. 

We ought to comment on two important things about 
the model: first of all, we should emphasize that, after 
spontaneous symmetry breaking, the Dirac mass terms, 
IlwiiiIr + vil^dvir, are not any more invariant under 
(if>r because of the non zero mass differences m T — 
and m T — m M , as we can see from Eq. (j3j|, whereas the 
Majorana mass terms do possess exact fj, f-> r symmetry. 
Second, the symmetry breaking in the diagonal Dirac 
mass matrix will generate the same effect in the effective 
neutrino mass matrix, M„, when the see-saw mechanic is 
implemented. Indeed, one gets the mass matrix 





( m ee 


m efl 






My = - 








(6) 




V m er 




m TT J 





that exhibits the breaking of the fj, •<-> r symmetry at tree 
level through the mass differences m eT — m e/J = yfh e (fh T — 
7^ and m TT — m w = w(m^ — m^) 7^ 0. 
Although, this symmetry is not exact in the neutrinos 
sector, it can be considered as an approximate symmetry 
as a result of a small breaking in the mass difference rh T — 
rh^. It is allowed, as we show later, because we have the 
freedom to manipulate the free parameters of the LRSM, 
specifically in Eq. (j3J). The effective mass terms can 
be separated as i> L (M v ) ^rU^ +Vl5M v v ( £, where the 
last term contains explicitly the breaking. Such a term 
affects lepton mixings, inducing calculable deviations in 
the mixing angles from the bimaximal mixing scenario. 



tan 20, 







sin 0i3 



sin a 



V8 



m ee ) 



1 


2m Mr <5 - 


■ em e/i 


V8 




- m MT - m ee ) 


I 


£ (™ W _ m »r - 


m ee ) + 2Sm efi 


~4 




- m^r - m ee ) 



(7) 

;(8) 

,(9) 



and e = m T 



0, 



where we have defined m eA1 
\(m TT +m w ), 5 = m eT ~m, 
last two parameters give us an account of how strong the 
breaking of the fi •<-> r symmetry is. We expect them to 



m w . The 



be small, compared to m ell and 



u , respectively, since 
they are related to a and $13 as we can see from Eqs. ^ 
and ([9]). On the other hand, the eigenvalues of M v are 
given as 



mi ? 


a m ee 


m 2 ? 




m 3 p 





V2 



(10) 

(ii) 

(12) 



Next, we would like to get some predictions out of this 
model, particularly for Eqs. {[sj) and ^ that depend 
on the elements of M v . In order to do this, we have to 
identify and properly choose the free parameters of our 
model. As we can see in M Vl if one writes this matrix 
and Eqs. ([7])-([9]) , in terms of Dirac and Majorana masses, 
one realizes that there are seven parameters: x, y, w, 
and z, which stand for the elements of the inverse matrix 
, and rh e pL T corresponding to Dirac masses. We 
would like now to fix some of those parameters using 
neutrino observables, others than those we are intending 
to predict at some level. Thus we take Am|, = m\ — m\, 



Am\ TM = \\m\ + m\ — 2m\ \ , the scale of neutrino mass 



hierarchy m 3 , and the solar mixing angle 9q as our inputs 
to fix the free parameters as much as possible. In spite of 
the lack of observables, we can only determine four of the 
model free parameters in terms of these known ones. So, 
we decided to fix x, y, w, and z, keeping in mind that the 
Dirac neutrino mass matrix is not invariant under ji •<-> r 
symmetry, and, as we have commented, that this is the 
only source of symmetry breaking. Thus, from the above 
given definitions of AttIq, Ato^ tm , and using Eqs. (|7j) 
and ( 12 ), we can solve the system of variables for the four 
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inverse Majorana masses, which will be now given as 

x = [sin 2 9 Q A + cos 2 9 Q B] ; 

sin 26» 
V2m e (m r + m^j 



where 



1 



(m 2 . + m 2 ) 



[sin 2 8 Q B + cos 2 6» Q A + m 3 ] 



2m T m M 



[sin 2 6» B + cos 2 6» A - m 3 ] , (13) 



.4 



m 2 ± Am^Af + -Am| ; 



m 2 ± Am| TM - -Am Q ; 



(14) 



between charged lepton and Dirac masses from Eq.Q. 
Thus, without loss of generality, we can take k t <C k% 
and assume that Yukawa couplings hi and h 2 are smaller 
than hi and h 2 , respectively. Then, we easily imply that 
the parameter h%jh 2 should be small to generate a tiny 
symmetry breaking. However, knowing that hi < hi, 
its ratio will be further constrained for the condition 
h s = b^/hx <§; 1, and as a result of the definition of r and 

R they will appear to be only dependant on the new pa- 
rameter h s - Therefore, with this simple approximations 
we can reduce the number of effective free parameters 
from seven to five. It is worth stressing that this is pos- 
sible since the final expressions for sin #13 and sin a do 
not depend on Dirac neutrino mass but throughout the 
mass ratios r and R, which can finally be simplified by 
hiding one of the parameters in an effective ratio, as we 
have seen. 



and + ( — ) stands for the inverted (normal) hierarchy 
case. These expressions would allow us to simplify our 
predictions for the mixings to write them in terms of 
neutrino observables, as we shall see next. So far, we 
still have three (Dirac mass) free parameters but, as it 
is easy to realize, the first Dirac mass, m e , shall not be 
involved in the determination of sin #13 and sin a, since it 
has a null contribution to M„ at the lower seesaw order. 
Therefore, as far as it matters for the mixings, we end 
up with two unknown parameters, which can be reduced 
to one by further considerations, as we will comment on 
next. Hence, by introducing the above expressions for 
the four Majorana masses into Eqs. ^ and ([9|, one gets 
the analytical forms 



{2r- R)(Ac 2 + B s 2 0) -TTi 3 (2r+i?) 
m 2 -m 3 (A + B) + AB 



*10{A-By -2R(Ac^ d + E 



){As'*0 + Bc J 'e- 



l^m 2 -m 3 (A + B) + AB^j 



(15) 

(16) 



where we have defined r = (m T — m ll )/(m T + m^) and 
R = (fh 2 — rh 2 ^)/{m 2 + m 2 ). Also, in these formulas, we 
have used a simplifying notation where s 2 9, c 2 9, and s29 
(s 2 29) stand for sin 2 6> , cos 2 9 Q and sin20 (sin 2 20 ), 
respectively. It is worth remarking how interesting a 
feature the parametric m e Dirac mass independence of 
M v becomes. Because of this fact, m ee , which fixes the 
strength of double beta decay, is very well determined by 
the observables, as one can notice from the first expres- 
sion in Eq. (13), which gives 



,A + cos 2 6» 0J B 



(17) 



that only depends on the unknown mass scale TO3 and 
the hierarchy. 

In addition, we have two (fh T and fh^) free parame- 
ters whose difference is the only source of fi r sym- 
metry breaking. Given the ratio forms r and R, we can 
introduce a further phenomcnologically motivated reduc- 
tion of the parameters, by observing that fh T — fh^ must 
be small to truly generate a tiny breaking of the sym- 
metry. In order to do that, let us take the relation 



IV. MODEL PREDICTIONS 

As a result of the approximations introduced above, 
the expressions for sin #13 and sin a given in Eqs. fll5]|16| ) 
are simplified into 



sin 13 



sin a 



m^h s m 



sin2(9 (A — B) 
2 -m 3 (A + B) + AB 
m 3 cos 2(9 (A — B) — AB 



m 3 {A + B) + AB 



(18) 



U9) 



where m = (m r — m^)/(m T + m^). We note that these 
expressions have a linear dependence on h s , and this will 
allow us to get a direct correlation between the angles 
that do not depend on h s , as we will show later on. Next, 
although we could obtain explicit expressions for each hi- 
erarchical case, we rather prefer to explore a more pic- 
torial representation of the parameter space, by showing 
the experimentally allowed regions for h s and the hierar- 
chy mass scale 1713, such that sin #13 and sin a are right 
below the current experimental limits. To be more ex- 
plicit, we consider the upper limit value for sin 6*13 to be 
of the level of la deviation, as given by the fits of cur- 
rent neutrino experiments data [51 [3] . We take the upper 
limit to be sin #13 < 0.16 for the sake of demonstration. 
Additionally, we observe that | sina| ~ || — sin 2 #23!, and 
thus sin a can be considered to be at most of the level of 
la deviation from the central value in sin 2 9 2 3, such that, 
at the end of the day we arc indirectly measuring sin #23- 
Here we assume the upper value |sina| < 0.07. Under 
these two premises, we can plot the corresponding sin0 13 
and sin a contour values in the m 3 — h s parameter space, 
respectively, that should bound the allowed parameter 
space. Our results are depicted in figures 1 and 2. 

From Fig. [l] for sin #13 all the region below the curves 
is physically acceptable, within lcr deviation. At the 
same time, one observes that the sin# 13 allowed region 
is mostly insensitive to hierarchy, meaning, we can not 
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FIG. 1: The la region for sin #13 in the parameter space of 
h s and 777,3 , for both normal (solid lines) and inverted (dashed 
lines) hierarchies. The upper and lower contours have a value 
of sin 0i3 = 0.16, within a a deviation. Notice that solid and 
dashed lines completely overlap 



FIG. 2: The la region for sin a as a function of h„ and 777,3, 
for both normal (solid line) and inverted (dashed line) hierar- 
chies. The upper and lower contours correspond to the limit 
value of sin a = 0.07. 



differentiate between normal and inverted hierarchy. In- 
deed, the regions basically overlap, so that, both hierar- 
chies share most of the same parameter region. In par- 
ticular, two extreme points we have to comment on: for 
small values of h s , 777,3 is close to its maximal value. On 
the other hand, for large values of h s , 7773 is larger than 
its middle maximal value. In addition, from the figure 
one may identify a window of acceptable values, which 
lies within h s « (0.04 - 0.1) and m 3 w (0.26 - 0.4). The 
important result here is that h s must be tiny in order 
to have a consistent value for sin# i3 , according to the 
experimental limits. The same conclusion will arise for 
sin a, as we will see from Fig. [2j 

Before starting the discussion about Fig. [2] we would 
like to comment that we have taken the absolute value 
of sin a just for the sake of depicting together all cases 
in the same plot. However, the case of inverted hierar- 
chy produces negative values for sin a, which already in- 
troduces a way to distinguish the hierarchy if one could 
identify in what direction #23 deviates from the maxi- 
mal value. We also observe that normal and inverted 
hierarchy can easily differentiate at small values of 7773. 
However, for tiny values of h s (and large values for 7773) 
both of them are difficult to point out from each other 
(apart from the overall sign). This plot is more con- 
strained since the allowed parameter region for sin a re- 
mains below the contours for each hierarchy. Again, 
from this figure we infer the phenomenologically valid 
interval of values that gets within h s s» (0.001 — 0.07) 
and 7773 s» (0.001 — 0.22) for inverted hierarchy, and 
h s fa (0.001 - 0.06) and m 3 « (0.06 - 0.22) for normal 
hierarchy. Comparing the parameter regions of 7773 for 
sin #13 and sin a graphics, this seems to be in disagree- 
ment since both sets of values should be common in order 



to satisfy simultaneously such expressions. Nonetheless, 
in Fig. 3 we show the valid region of values for 7713 to 
both hierarchy cases. 

One interesting prediction of the model is that the ratio 
sin a/ sin 0i3 does not depend on the parameter h s , as 
we already pointed out, and thus, an appealing formula 
arises which represents the key expression that allows one 
to falsify this model 



sin a 
sin6>i 3 



cot 29, 



AB 







7773 sm29 Q {A-B) 



(20) 



This ratio is principally driven by the second term which 
has a clear dependence on 7773 and the rest of the observ- 
ables, and also an extra contribution comes from the so- 
lar angle. This can be enhanced or decreased depending 
on the hierarchy scheme, in particular, for strict normal 
hierarchy one obtains a fit value of about 68.45 when cen- 
tral values are taken for the observables. On the other 
hand, inverted hierarchy has a notable dependence on 
7773 whose value is not well determined. One can in gen- 
eral depict the expected ratio in terms of 7773, as we show 
explicitly in Fig. [3] 

As we can see from Fig. |3j there is a clear difference 
between both hierarchies for small values of the 777.3 mass 
scale. For normal hierarchy the ratio decreases with the 
hierarchy scale 7773 , whereas for inverted hierarchy the ra- 
tio rather increases (negatively) with 7773. However, it is 
hard to differentiate between the two hierarchies for large 
7773. The most important point we do want to make about 
this plot is that our model, under the previous approx- 
imations, predicts that due to the correlation between 
I sin a I and sin 13 , the latter parameter will very likely lie 
in below a range of 10 -3 , which can be difficult to reach 
in future experiments, given that | sina| ~ 1 5 ~ sin 2 #23! 
is already bounded below 10 -2 by current experimental 
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FIG. 3: The la region for | sin a>\/ sin #13 as function of 7713, 
for both normal ( solid lines) and inverted ( dashed lines ) 
hierarchies. 



data. Nevertheless, this correlation will make it possible 
to falsify our model in the following way. If future exper- 
iments achieve to measure |sina| and/or sin #13, then, 
they should lie in the allowed region according to Fig. [3] 
On the contrary, if such were not the case, the model 
would be automatically discarded. 



V. CP PHASE CONTRIBUTION 

Next, let us add some comments on CP phase contribu- 
tions. As we already stated, and as we will discuss below, 
it seems unlikely that such contributions may drastically 
change our results. In particular, aside from some small 
numerical changes, it does appear yet difficult that CP 
phases may increase the tight upper bound the model in- 
dicates for sin 6*i3, which has been put on at about 10~ 3 , 
and which is already far away from the current experi- 
mental bound (~ 10 _1 ). 

In order to show our present point in a simple way, let 
us briefly reconsider our analysis of Sec. III. As usual, 
the mixing matrix is in general given by U = UpmnsK \ 



diag{\, 



oi<Pl f>i4>: 



2 ) stand for the Majorana 



where K 

CP phases. Also, Upmns contains (in the standard 
parametrization) the ip Dirac CP phase. So that, under 
the same considerations as for Sec. Ill, M v is diagonal- 
ized by U, and Eqs.(8|9) get replaced by 



1 

7^ 



P* (S-eC)-K(S-eCy 



1 



4m„ 



P 2 - K 



' ( ft 



(21) 
(22) 



where, as before, 
the short-hand 

K : 



sinf9i3e l(p , and we have introduced 



notation 
and C = 



where P = 
= m eM /2m MT . 



"'/it- «>">> >..• — "'.///-' "7 - r- Here, 5 and e pa- 
rameters have the same form as in previous sections. We 
must keep in mind the M u matrix elements are now com- 
plex (| m a p I e lPj , where a,/3 — e,fi,T, and j — 1,..,6), 
and the Majorana CP phases are implicitly included in 
Eqs. ( |21p2l ). Of course, when CP phases are canceled 



(taking all masses to be real), the above expressions re- 
duce to the previous results as given in Eqs. (|8j9]). In- 
teresting enough, Eq. ([7]) maintains its form even in the 
presence of CP phases. 

Having written the modified formulas for sin #13 and 
sin a, we are now interested in figuring out how CP- 
violating phases affect our previous considerations and 
conclusions. Performing a complete and general analysis 
of the CP-violating case is challenging and complicated. 
Nevertheless, there is a simpler way to proceed for our 
purpose, and that is to use the numerical results that 
came from our previous analysis as a grounding point, to 
estimate, through a general example, the effect of reintro- 
ducing non-trivial CP phases. Hence, we will proceed as 
follows: first, we use our CP-conserving results to rebuild 
the effective neutrino mass matrix elements, thus identi- 
fied as |TO a /3|, for a given set of initial inputs, and then 
reintroducing CP phases on M„, and recalculating our 
predictions to determine how far the new results are from 
those expected from the CP-conserving case. We will 
particularly focus on the correlation given in Eq. ( |20[ ), 
which is our main result here. Moreover, since the abso- 
lute values of bjra etl and e/m MM are supposed to be small, 
we may conclude that m e ^ and m er will have about the 
same phase, as so will to mm and m TT , irrespective of 5 and 
e own arbitrary phases. For simplicity, we will suppose 



that 



and S have the same phase, and analogously 



for m^fj, and e; this choice shall reduce the number of 
CP phases on M„ to four, although, it is clear not all 
will be truly physical. Under these considerations, the 
correlation between sin a and sin #13 is now expressed as 



sin a 



1 



\z\ ^/2\m^ T m ee \ 
where we have defined 



h 



(23) 



fi 



1 1 — 2 1 



o l 1 



Wl 



and 



h = 



0' 



(24) 



(25) 



Here, 7 = pi+ p 3 — 2p 2 and 71 = p\ — pz stand for the two 
effective CP-violating phases entering in our results that 
came out of combinations of the arbitrary phases we first 
introduced in the neutrino mass matrix. Therefore, we 
realize the correlation among mixings depends only on 
these two CP phases, apart from the 7713 hierarchy mass 
scale. 

Next, we can use Eq. (23) to quantify the phase de- 
pendence of our previous CP invariant results. This we 
do by exploring some general examples, by taking spe- 
cific values for 7 and 71, and then replotting the ratio 
|sina|/|z| as functions of m 3 for the central values of 
the observables. We show in Fig [4] (solid line) the re- 
sults we obtained for the CP invariant case, that is for 



7 = and 71 = 0, for both neutrino mass hierarchies, 
as before (see figure [3] for a comparison). Next, in same 
figure, short-dashed lines have been obtained when the 
arbitrary values 7 = 7r/4 and 71 — it/ 5 arc taken. Sim- 
ilarly, for 7 = it/ 2 and 71 — tt/A the associated plot is 
represented by the dotted line. Finally, the long-dashed 
line stands for the arbitrary set of values 7 = ir/7 and 
7i = tt/2. 

\ssBm\ 
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FIG. 4: Upper bound for |sina|/|z| as a function of 7723, for 
different values of 7 and 71 as discussed in the text. As before, 
for normal hierarchy, the contours (the correlation) decrease 
for small values of 7713 whereas for inverted hierarchy the con- 
tours increase (negatively). 



From Fig. Q, one may deduce important facts about 
the CP- violating phases contributions to the correlation 
between sin a and sin 6* 13 : In general, for both hierarchies, 
CP phases do not change the general profile of the la 
regions, and neither seem to introduce drastic changes 
on the sin aj sin #13 ratio, which at least for the arbitrary 
values we have explored, keeps sin #13 upper bound on 
about 10 -3 , as we have pointed out before. Indeed, in 
spite of considering the largest changes we found on the 
mixing correlation (for 7 = tt/7 and 71 — tt/2), the result 
is not yet enough to significantly reduce the bound. 

Nonetheless, as a cautionary word, as we have not pre- 
sented here a truly general analysis for the CP-violating 
case, we can not discard the possibility that CP phases 
may conspire among themselves to produce an accidental 
enhancement of the sin a/ sin #13 ratio. A naive analysis 
shows that this seems to happen for instance for the pre- 
cise values of 7 — tt and 71 — it, which set the mentioned 
ratio to the order of 10 , which in turn would locate 
sin #13 close to the experimental sensitivity level. How- 
ever, such a case has almost no parameter space on the 
CP phase sector, and it seems very unlikely that such 
extreme values could be potentially physical, since even 
phases are subject to higher order corrections, and no 
symmetry exists that may protect phases as to keep them 
at those precise values. Thus, a conservative point of 
view would point towards the strength of our previous 
CP invariant results. A more careful and general study 
of CP violation in/if>T symmetric models is in progress, 
and a detailed analysis of this scenario will be presented 
soon elsewhere. 



VI. CONCLUSIONS 

We studied the lepton sector of the LRSM dressed with 
a (Z2) 3 discrete symmetry, as was done in two frame- 
works: (1) neglecting the CP- violating phases and (2) 
involving them. Using discrete symmetries we build a 
framework where the charged lepton mass matrix is di- 
agonal at tree level, so that we start since the beginning 
in the true flavor basis, and get a model that naturally 
realizes a slightly broken fi t symmetry. In addition, 
we have determined #13 and the deviation from the max- 
imal value of #23 in the limit of a small fj, •<->■ r symmetry 
breaking that appears at tree level in the effective neu- 
trino mass matrix. These two expressions are fixed in 
terms of one free parameter of the model, h s , the mass 
scale hierarchy 777,3 (for the first scenario), and two CP- 
violating phases (for the second one); within the first 
scenario, a set of values for h s and 7713 were found which 
are consistent with current limits imposed by the exper- 
imental data for sin #23 and sin #13 . 

The present model is quite elaborate, but predictive. 
In particular, a correlation among sin a that measures the 
deviation of sin 2 #23 from 1/2, and sin #13 arises which 
only depends on one neutrino mass scale, 7713 (and the 
CP phases). This is a remarkable result and provides 
a way to test the model, since the observed value of 
I sina|/ sin #13 would also provide a direct determination 
of the scale of 771,3 for a given hierarchy. We should point 
out this ratio does predict an allowed region for sin #13 
which must be below 1CU 3 , in order to satisfy the cur- 
rent experimental data, that upper bounds the expres- 
sion 1/2 — sin 2 #23- Such a sensitivity could be hard to 
reach in future experiments but it seems possible. On 
the other hand, as an extra bonus of this model, a clear 
determination of the m ee strength of double beta decay 
was found. 
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VII. APPENDIX 

The most general form of the scalar potential of the 
LRSM was analyzed in [9], and it has 18 terms which 
are invariant under the gauge group SU (2) L x SU{2) R x 
U(1)b-l and parity. Now, from our scalar particle con- 
tent: two bidoublets <pi (i — 1,2), and four triplets, 
Aj (j = 1,...,4), we should write the potential. Even 
though, it is larger than the original potential, it is pos- 
sible to show by a lengthy calculation that it has a least 
a minimum. Here we are interested in showing that the 
triplets A2R and A3A have the same VEV when we mini- 
mize the potential. So, without loss of generality, we just 
write the terms in which their VEV'-s are involved after 
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spontaneous symmetry breaking and they are given as: 



V 



n\ + 2bv 2 R + 2b 1 v\ R + ibakik! 



+b 2 (kl + k{ \ + + h \ k\ + &2 2 ) + 46 6 fc 2 fc2 



-b 7 k' 2 2 



b 8 (v. 



2R T v 3R J 



2bqV, 



2 v 2 



9 V 2R V 3R 



(26) 



We can easily conclude from the above expression that 
A 2 r and A 3 r have the same VEV, by directly calculat- 



0. 



ing the two minimization conditions, £ v — £ v 

o _ dv 2 R ov 2 r 

It is not difficult to realize that this is a consequence of 
the fact that V remains symmetric under the v 2R <H> v 3R 
exchange, derived from the original Z 2 symmetry. This 
ensures that the Majorana mass matrix, Mr, will remain 
invariant under fi -H- r symmetry after spontaneous sym- 
metry breaking. On the other hand, we have to keep 
in mind that all VEV'-s will only depend on potential 
parameters after solving the simultaneous conditions for 
each VEV. 
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